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Abstract 

A weak law of large numbers is established for a sequence of systems of N classical point particles 
with logarithmic pair potential in R n , or S™, n € N, which are distributed according to the 
configurational microcanonical measure 5{E — H), or rather some regularization thereof, where 
H is the configurational Hamiltonian and E the configurational energy. When N — > oo with 
non-extensive energy scaling E — N 2 e, the particle positions become i.i.d. according to a self- 
consistent Boltzmann distribution, respectively a superposition of such distributions. The self- 
consistency condition in n dimensions is some nonlinear elliptic PDE of order n (pseudo-PDE if n 
is odd) with an exponential nonlinearity. When n = 2, this PDE is known in statistical mechanics 
as Poisson-Boltzmann equation, with applications to point vortices, 2D Coulomb and magnetized 
plasmas and gravitational systems. It is then also known in conformal differential geometry, 
where it is the central equation in Nirenberg's problem of prescribed Gaussian curvature. For 
constant Gauss curvature it becomes Liouville's equation, which also appears in two-dimensional 
so-called quantum Liouville gravity. The PDE for n = 4 is Paneitz' equation, and while it 
is not known in statistical mechanics, it originated from a study of the conformal invariance 
of Maxwell's electromagnetism and has made its appearance in some recent model of four- 
dimensional quantum gravity. In differential geometry, the Paneitz equation and its higher 
order n generalizations have applications in the conformal geometry of n-manifolds, but no 
physical applications yet for general n. Interestingly, though, all the Paneitz equations have an 
interpretation in terms of statistical mechanics. 
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1 INTRODUCTION 



Early in his scientific career Joel Lebowitz spent a postdoctoral year or so working with Lars 
Onsager. There are many scientific similarities between these two towering figures of statistical 
physics: both in possession of penetrating mathematical powers; both broadly interested in physics; 
both with a good taste for choosing problems worthy to work on; and both with deep insights into 
physics that have inspired - and continue to inspire - many other physicists and mathematicians 
around the world. Of course, Joel has also worked on many problems that were inspired by some 
work of Onsager, and while it would be impossible for me here even to merely mention them all, one 
contribution to a problem of Onsager that I was privileged to collaborate on with Joel has recently 
had some interesting mathematical spin-off, and I am very glad to present these new results in this 
paper in honor of Joel's 70 th birthday. 

In the next section, I will first describe Onsager's application of statistical mechanics to point 
vortices in two dimensions, the questions it has raised, and our contribution to it. In the two 
subsequent sections I will describe Nirenberg's problem in the conformal differential geometry of 
two-dimensional manifolds and its recent generalization to n-manifolds. Point vortices in two 
dimensional flows have logarithmic pair interactions. In the last section I will explain how statis- 
tical mechanics of logarithmic interactions, but now in all dimensions, provides solutions to the 
Nirenberg-type problems in conformal geometry. 



ONSAGER'S POINT VORTEX DISTRIBUTIONS 



In a pioneering paper on statistical fluid dynamics |57|| , Onsager studied, among other things, 
the microcanonical ensemble of N classical point vortices in a two-dimensional, incompressible 
Euler flow. The motion of iV vortices in a connected domain A C R 2 is governed by Kirchhoff's 



Hamiltonian [H3l 44 



^)(x 1 ,...,x iV ) = i^ 5>c,G*(x i ,x J -)+ ]T CiFWfr), (1) 

l<i,j<N l<i<N 

where Xj G A is the position of the i-th vortex; q its circulation, expressed as dimensionless multiple 
of a suitable reference unit; G*(x, y) = G*(y,x) is a renormalized Green's function for —A on A, 
i.e. G*(x, y) = G(x, y) if x ^ y, with G(x, y) solving — AG(x, y) = 27r£ y (x) in A, and 



G* (x, x) = Jim (G(x, y) + In |x - y | J . (2) 

In case of a bounded domain with piecewise regular boundary dA, the physically appropriate 0- 
Dirichlet boundary conditions for G are imposed on dA, but in principle other boundary conditions 
are of interest too. Finally, F^ N ' is some externally applied stream function (whose strength may 
be proportional to N). Considering a bounded domain with finite area |A|, and noting that up to 
a trivial factor v/jc^] the canonically conjugate variables of the vortices are given by the Cartesian 
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components of their positions in A, Onsager observed that the phase space volume of the set 
jijW < E} in the ./V vortices phase space A^, 

$A(E,N) = \{hW <E}\ , (3) 

is a monotonically increasing function of E, bounded by |A|^. He noticed that this implies that 
Boltzmann's entropy S A (E,N) = In (\A\~ N <fr' A (E, N)) (where & A (E,N) = d E <S>A(E,N)) must 
reach a maximum at a particular value E m (N;A) of energy, such that |5^] , p. 281: "negative 
"temperatures" ...will occur if E > E m , ... vortices of the same sign will tend to cluster ...[and] 
large compound vortices [be] formed in this manner..." 

Onsager's insight not only predated Ramsey's prediction of negative temperatures in spin sys- 
tems [33 1 by 7 years, Onsager was (once again) way ahead of his time. His program was not picked 
up until almost a quarter century later, when Montgomery pointed out 5C] that numerical sim- 
ulations of turbulent fluid flows suggested an explanation in terms of Onsager's 1949 paper. But 
Onsager had not made any attempt to extract continuum vorticities from his statistical vortex 
distributions, and so several authors [36, 51, 61, 46, [3|] now came up with the proposal to use 
mean-field theory for this purpose. 

As it is with traditions, and statistical mechanics surely has a long tradition, certain general 
wisdoms tend to be passed on which sometimes may not be so generally valid. One such general 
wisdom holds that 'mean- field theory is wrong.' This traces back to the failure of the prototype 
mean-field theories of van der Waals and Curie- Weiss to give the correct data for the critical 
point of the condensation and magnetization phase transitions, respectively. In particular, they 
predict incorrect critical exponents. However, as I recall Michael Fisher pointing out in a lecture 
on Coulombic criticality, mean-field theory is not that bad after all. Another general wisdom states 
that thermodynamic concepts such as temperature become precise only in the thermodynamic 
(bulk) limit of an infinitely big system. A third general wisdom says that the statistical ensembles 
are equivalent. Against this background it is easier to appreciate that it took a while to realize 
that neither wisdom is correct in the case of Onsager's statistical theory of vortex clustering. 

Indeed, the first attempt [^lj to put Onsager's theory on a more rigorous basis was made in 
terms of the traditional bulk thermodynamic limit for a neutral two-species vortex system, using 
formal central limit arguments of Khinchin which where asserted to apply also to the negative 
temperature regime. Moreover, in that paper [51| also the BBGKY hierarchy was considered, and 
the mean-field equation for the distributions at negative T, which in []36| had been obtained by the 
standard van-der-Waals-theory type combinatorics, was now obtained as Vlasov approximation. A 
subsequent rigorous study by Frohlich and Ruelle [33] however showed that negative temperatures 
do, in fact, not exist in the bulk thermodynamic limit of a neutral two-species vortex system. 

Inspection of the phenomenon that Onsager predicted reveals, however, that O(N^) truly long 
range pair interactions in the k th cluster of vortices of the same sign are involved |46]| . Hence, we 
are dealing with a strictly nonuniform system with non-extensive energy scaling. Understood from 
this perspective, the relevant limit N — > oo in which Onsager's prediction of negative vortex tem- 
peratures attains a rigorous meaning is not the bulk thermodynamic limit but rather a continuum 
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(Euler) fluid limit. Moreover, in this limit mean-field theory should become exact in the sense of a 
weak law of large numbers. 

Thus, the following picture emerges. For the neutral two-species system, with \a\ = 1, and with 
F = for simplicity, distributed in a bounded domain A C I 2 according to the microcanonical 
measure 

^ N ' E \d Xl ...dx N ) = ^ ^y «5 (E - H^) d Xl ...dx N , (4) 

where dx denotes Lebesgue measure on R 2 , we have to fix A C M 2 and e = E/N 2 > 0. Then, in 
the limit N —* oo, the Boltzmann entropy per vortex will converge to a continuous function of e, 

l^±ln(\A\- N $' A (N 2 e,N))=s A (e), (5) 

which is given by a variational principle, 

s A (e) = max £ - / p° \n(\A\p°)dx , (6) 

P ,P a= ± J A 

where the maximization is carried out over the probability densities p^ which satisfy the energy 
constraint 

J G(x,y)u(x)u(y)dxdy = e, (7) 

where u> = p + — p~ . Moreover, the microcanonical equilibrium measure itself will converge, for 
e > and in a suitable topology, to a convex linear superposition of infinite products of those 
absolutely continuous one- vortex measures whose densities pf are the maximizers for ©-(0). In 
particular, if the superposition measure is a singleton (a unique maximizing pair pf exists), then 
we have the weak law of large numbers, that for all / S C°(A), 

N/2 

t)pt(*)dx (8) 



.7=1 JA 



in probability. In (||), the summation extends over either all positive or all negative vortices. If the 
superposition is not a singleton, pf in (||) is to be replaced by the corresponding superposition of 
probabilities. Furthermore, the maximizers satisfy the self-consistency conditions 



exp ( 



/ G(x,y)u; e (y)dy 

JA 



(9) 



where uj e = pf — p~ , and where the constants (3 and p^r h are to be chosen so that the constraint 
Ja P^dx = 1 and (0) are satisfied. Not only are the equations @ clearly mean-field, solutions having 
(3 < are known to exist, and the corresponding continuum vorticities oj £ satisfy the stationary 



Euler equations for incompressible flows [47|. Onsager's prediction of negative temperatures in 
vortex systems indeed attains a precise meaning in this limit. 

There can be hardly any doubt that the above picture is correct, and that it will eventually be 
verified rigorously. An almost complete verification is available for the mildly simpler single species 
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system. By generalizing a result of Messer and Spohn jSJ for Lipschitz continuous interactions to 
the logarithmically singular point vortex interactions, first the canonical point vortex ensemble was 
conquered, independently in || and [39|. While the microcanonical ensemble in the strict sense 
described above has not yet been treated, rigorous works exist in which the microcanonical point 
vortex measure 5(H — E) is replaced by a regularized measure, the regularization being removed 
after the limit N —* oo has been taken. Again, the first result, by Eyink and Spohn [pi]] , was 
for regularized interactions. The singular point vortex interactions where then treated by Caglioti 
et al. [ 10 1 . The limit N — ► oo in |1C] is, however, constructed under the assumption that the 
microcanonical and canonical ensembles are equivalent, an assumption which is not generally valid 
for these finite-domain mean-field limits. This last barrier was finally overcome by Joel and myself 
in [41], where we go beyond equivalence of ensembles for the logarithmically sing ular §). Of course, 
the mean-field theory obtained is precisely the single-species version of the one described above. 

Future work should extend these results to the neutral two-species vortex system. In contrast 
to the many rigorous results obtained for these neutral systems in the traditional thermodynamic 
limit, see |33| and references therein, only few facts are rigorously known for their Euler fluid 
limit. Curiously enough, while Frohlich and Ruelle |3^[ proved absence of negative temperatures in 
the standard bulk limit, in |4^[ I was able to prove absence of positive temperatures in the Euler 
fluid limit for the neutral two-species system; in fact, (3 < is bounded away from zero. This 
raises the question as to the whereabouts of Onsager's E m , at which temperatures are supposed 
to switch from positive to negative values in a neutral system? The answer is that there is much 
room between the low-energy regime E = Ne of |33|| and the high-energy regime E = N 2 e of 
pT| , [U], |31| . In fact, in [11] it was found that for a neutral two-species system E m (N) ~ CNlnN, 
falling precisely inbetween these two regimes. The vicinity of E m , analyzed in |55], O, turns out 
to be a small-entropy regime where S, not S/N, converges to a limit when N — > oo. 

The results in [55, 11] are the only ones obtained directly from 5{H — E). The construction of 
the Euler fluid limit N — » oo directly from 5(H — E) remains an open problem. For now, [41] is the 
final word in the construction of the Euler fluid limit for Onsager's statistical vortex distributions. 

Interestingly, 41] also holds the key to some answers to questions in conformal geometry, to 
which we turn next. 



3 NIRENBERG'S PROBLEM 

Many years ago, see |^3|, Louis Nirenberg raised the following question: "Which real functions K(x) 
on § 2 are Gauss curvature functions for a surface over § 2 whose metric g is pointwise conformal 
to the standard metric go on § 2 ?" To see what this question has to do with Onsager's vortex 
distributions, we need to recall a few basic concepts of differential geometry as found, for instance, 
in |66[| . As we shall see, the connection with statistical mechanics, or at least thermodynamics, 
must have been suspected by differential geometers long ago! 

Consider first, for simplicity, surfaces over M 2 that are embedded in M 3 . Thus, let x E A C M 2 . 
The two Cartesian coordinates x l ,x 2 of x provide two independent real parameters for the local 
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represention of a surface S C M 3 , given by S = X(x), x G A, with X G M 3 . The line element 
on S is given by da 2 = gijdx % dx^ , where we use Einstein's summation convention, and where 
9ij = {^x i -X,d x jX) (Euclidean inner product) are the components of the metric tensor. If v is the 
unit normal at S induced by the representation X(x), then the Gauss curvature K(x) is defined 
by (x means cross product in M 3 ) 



d x w x d x 2V = K{x) d x iX x d x iX 

Gauss' Theorema Egregium asserts that K(x) depends only on the gij. Writing gn 
and (722 = G, and moreover x 1 = s and x 2 = t, we have the Frobenius formula 



K 



+ 



(10) 

E, 512 = F, 




(11) 



see |3|, 56|, from which it follows that K is independent of the parametric representation X(x) of S. 
This freedom can be used to simplify (pd|). In particular, the representation for which E = G and 
F = is conformal, i.e. da 2 = Edx 2 . Interestingly, it is known in differential geometry as isothermic 
parameter representation. In this representation, (11) reduces to K = —-^ [d 2 \a.E + df hiE), or, 
setting InE = 2u and recalling (s,t) = x, 



K{x) = -e- 2u(x) A«(x) 



(12) 



which gives us K when u (i.e. E) is given. Nirenberg's question, here its analog on M 2 , addresses the 
inverse problem, i.e. to prescribe the putative Gauss curvature function K and study fll2]) as a non- 
linear elliptic PDE for the unknown u(x), x G A C M 2 . In particular, when K is constant, K = ±1 
by scaling, then (|i~2| ) is known in differential geometry as Liouville's equation p5[| . Clearly, Liou- 
ville's equation is identical to what in statistical physics goes under the name Poisson-Boltzmann 
equation, 

- A^(x) = 2Tie^ h -^ ] (13) 

for the spatial density of a 'perfect gas' in A C M 2 , in thermal equilibrium at temperature f3~ l and 
chemical potential [i c h, distributed in its own Coulomb (or Newton) potential field tp. Hence, the 
notion of 'isothermic parameters,' so it seems; however, I have not been able to trace the originator 
of this terminology. In any event, (|l3|) is of course identical to the one-species specialization of @, 
i.e. (||) with p + = p and p~ = 0, and with V'(x) = J\ G(x, y)p(y)dy the stream function, so that 
— At/j(x) = 27rp(x). We have come back full circle to Onsager's problem. 

Liouville p5| himself already showed that the equation derived by him is in a certain sense 
completely integrable. In particular, all entire solutions on R 2 with finite integral curvature k = 
J" K 2 Ke 2u ^dx have been identified; namely, no entire solution exists when K = — 1 58], while 
for K = +1 Liouville's equation is solved when exp(2u) is any point on the conformal orbit of the 
Jacobian of the stereographic map S 2 



£ 2 , and k = 47r for all solutions, then; see p3| , 12, ^7], 21] 
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For the prescribed Gauss curvature problem in all M 2 with K ^ const., no general solution of 
( |j~2|) is available, yet over the years a vast knowledge has accumulated, e.g. |5f|, [54], 48], ||], [25|, 
[24], [20], [26], [22], |62|]. Of particular interest are radially symmetric Gauss curvature functions, 
for then the interesting question arises whether solutions u exist that break the radial symmetry. 
In the following, the sign of K (as a function) is defined as: sign (K) = 1 if K ^ and K(x) > 
for all x; sign (K) = — 1 if K ^ and -fT(x) < for all x; and sign (K) = if K = 0. In all other 
cases sign (K) is not defined. The next theorem is taken from |22j . 

Theorem 1: Assume K £ C Q (R 2 ) is radially symmetric, has well-defined sign, and satisfies 



[ \K(x)\e 2h(x) \x\ x dx <oo 



for some non-constant, harmonic function h : 



and some A > 0. Assume also that 



[ \y- x\-"i\K(x)\e 2H( - x Ux — >0 as|y|^oo, 



(14) 



(15) 



for all < 7 < 2. //sign (K) = —1, define 

K *(K;h) = -2?rsup{A > : (14) is satisfied} (16) 
where k* might be —oo for some K < 0. Then, for any such K and h, and k satisfying 



k e < 



' (k*,0) if sign (K) = -1 

[0] if sign {K) = 

(0, 4vr) if sign (K) = +1 



(17) 



there exists a classical solution u = Uh^ K of (\lq) with prescribed Gaussian curvature function K , 
uniquely if K < 0, having integral curvature 



K{yi)e 2Uh >^dx 



and asymptotic behavior 



Uh K (x) = h(x) - — In |x| + o(| In |x| 

Z7T 



as x 



oo . 



(18) 



(19) 



Theorem 1 considerably generalizes an earlier result of this kind (Thm.III in |25|]) which is restricted 

to compactly supported K in M 2 . Theorem 1 is proven in 22] as corollary of the construction, also 
given in |22[ ], of the Euler fluid limit in R 2 for the canonical point vortex ensemble, which generalizes 
the one for finite domains AcM 2 J39| , [9[. The canonical point vortex measure on M. 2N is given by 

1 



(N ' (3) {dx 1 ...dx N ) 



Z(J3,N) 



exp 



-pN^H^ dxi...dx N , 



(20) 



with G(x,y) = -ln|x-y| and fW(x) = -NP'^ln \K (x)| +2^t(x)] in H^ N \ (|). The reciprocal 
temperature (3 in ( |20| ) and the integral curvature k in Theorem 1 are related by k = —(3ir, with 

-/?€ («7M)- 
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Encouraged by these achievements of the canonical statistical mechanics approach to the pre- 
scribed Gauss curvature problem on M 2 , we now return to Nirenberg's problem in its original setting 
on the sphere E> 2 . Writing the conformal deformation of the metric as g = e 2u ^ go , where go is the 
standard metric on 8 2 , the problem is to find all K(x), x G § 2 , for which 

- Au(x) = K(x)e 2u{x) - 1 (21) 

has a solution u on S 2 . Here, A is the Laplace-Beltrami operator on S 2 w.r.t. go. Superficially 
(21) is hardly any different from the prescribed Gauss curvature equation on R 2 (|T2]). However, 



Nirenberg's problem for S is a hard problem, indeed; see [53|, |37 ], |15 |, |l6 l, [pjjl ; see also [37| 



[29] for more general compact Riemann surfaces. 

There are many obstructions to finding admissible K for (|2l|). The celebrated Kazdan- Warner 
theorem [[37j is one of them. According to |57j] a Gauss curvature function K on §> 2 cannot be 
axially symmetic and monotonic, unless it is a constant function, say K = 1. In the latter case 



the problem is completely solved, 52], jy], 56]. Another, more serious obstruction is the famous 
Gauss-Bonnet theorem [w| , which relates the Gauss curvature K g on a general compact 2-manifold 
(M 2 ,g) without boundary to the Euler-Poincare characteristic x(^ 2 ) 5 a topological invariant, by 

X(M 2 ) = ^ / KgdVOlg . (22) 



M 2 



Since x(§ 2 ) > 0] it follows from (22) that K(x) has to be positive for some x G S 2 , [ 53 1 . Furthermore, 
we have x(§ 2 ) = 2- Thus, k = An, i.e. 

' K(x)e 2u(x) dx = 4Tr, (23) 



which of course results also by directly integrating (^) over § 2 Contrast (||) with the range of 
integral curvatures covered in Theorem l! In particular, notice that k = 4-7T is not included in 
Theorem 1, and this is a matter of principle. Indeed, the restriction to —/3 < 4 in the canonical 
ensemble is an obstruction imposed on us by the local integrability requirement of the singularities 
in (|2"0|), and this will not improve if we put our vortex system on the sphere [42]. Thus, (p3) dashes 



any hope to apply the Euler fluid limit of the canonical point vortex ensemble at j3 = —4 to the 
prescribed Gauss curvature problem on S 2 . 

Not all hopes are dashed, though. Physically speaking, what happens at —(3 = 4 in the finite N, 
single-species canonical ensemble is that the vortex system concentrates onto a single point. This 
is done at the expenses of the 'heat bath,' which delivers a positive infinite amount of energy into 
the vortex system at fixed negative temperature. Taking the limit iV — > oo of this concentrated 
singular state gives a Dirac 5 mass on § 2 , corresponding to e = +oo. However, performing the 
mean-field limit iV — > oo first for —(3 < 4 produces a regular solution of the resulting analog on 
S 2 of the Poisson-Boltzmann equation (|l3|), and subsequently letting —(3 /* 4 may, or may not 
result in a regular limiting solution at —(3 = 4. This can be made more precise with the so-called 
concentration-compactness alternative of P.L. Lions, compactness of the sequence of solutions as 
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—(3 y 4 implying a regular limiting solution at — (3 = 4. The canonical ensemble may never- 
theless not provide enough control to decide the concentration-compactness alternative. However, 
whenever compactness holds, a regular solution exists which has finite energy e, and some finite 
entropy s. We conclude that the study of the mean- field limit for the microcanonical ensemble of 
point vortices on § 2 , with external stream function F^ N \ can be expected to yield valuable new 
information on Nirenberg's problem. 

As remarked earlier, the microcanonical ensemble in its strict sense has not yet succumbed to 
rigorous analysis. However, for the prescribed Gauss curvature problem on § 2 the construction 



given in [41] is fully sufficient. As a matter of fact, we shall state our microcanonical results in the 
more general setting on § n , n > 2. To see what this now is about, we need to briefly digress into 
Paneitz theory. 

4 PANEITZ EQUATIONS 

Initiated by the conformal covariance in Minkowski space-time HL 1 ' 3 of the Maxwell equations of 



electromagnetism, S. Paneitz [59| in 1983 discovered a quartic conformally covariant differential op- 
erator for arbitrary pseudo-Riemannian 1, 3-manifolds, together with an associated new conformal 
curvature invariant. Recently this theory has received considerable attention, see ||, ||, |fT?j] , @, 



[13], [14]. The significance of the Paneitz curvature, Q g , becomes apparent through an analog of 
the Gauss-Bonnet formula for a general compact 4- manifold (M 4 , g) without boundary, |l£ 



*( m4 ) = ^2 f M4 Ql^(*)| 2 + Qs(x)) d™l g (24) 

Here, x(-^ 4 ) is the Euler-Poincare characteristic of M 4 , and W is the pointwise conformally in- 
variant Weyl tensor. Moreover, Q g is defined by 

6Q g (x) = -A g R g {x) + ~R 2 g (x) - 3|Ric~ 9 (x)| 2 (25) 

where A g is the Laplace-Beltrami operator on (M 4 ,g), R g is the scalar curvature and Ric ff the 
traceless Ricci tensor, see p9|, |l3|]. Associated with Q g is Paneitz' quartic conformally covariant 
operator, such that, given (M 4 ,go) an d a conformal change of metric written as g = e 2u ^go, 
x G (M , go), the new curvature Q g is given by 

Q g (x) = e~ 4u ^ (A 2 u(*) + 5 go (^R go (x)I - 2Ric 90 (x)) d go u(x) + Q^x)) (26) 

Here, d is the differential and 5 the divergence. 

An analog of Nirenberg's problem on the 4-manifold (§ 4 , go); with go the standard metric, can 
be formulated thus: "Which real functions Q(x), x G S 4 , are Paneitz curvature functions for a 
4-manifold whose metric is pointwise conformal to the standard one?" This analog of Nirenberg's 
problem can be rephrased in terms of d26|), namely to find all functions Q(x) on S 4 such that the 
fourth-order PDE 

A 2 u(x) - 2Au{x) = Q{x)e iu{ - x) - 6 (27) 
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has a solution u on S 4 . Generalizations to S n (and, hence, to W 1 ) of the Paneitz equation (27) have 
been derived as well. On S n , n > 2, we have 



P n u{x) = Q{x)e nu ^ - (n - 1)! 



(28) 



with 



Pn= < 



re — 2 
2 



j^J (—A + A;(n — k — 1)); n even 



fc=o 



-A + 



n — 1 



\ 2 ~ 

J J| (—A + k(n — k — 1)); n odd, 

' k=0 



(29) 



sec @, ||, ||], 01 . On M n , we simply have 

(-A) n / 2 u(x) = Q{x)e nu{x ^ 



(30) 



with x G M n now. The operator P n in ( |28|) is the Paneitzian, and Q(x) the Queervature (pardon 
the pun) of order n. Notice that for n odd, P n is a pseudo differential operator. 



5 LOGARITHMIC INTERACTIONS IN ALL DIMENSIONS 



The increased complexity of the operators P n given in (|29| ) for high dimensions gives the Paneitz 
equations ( |28| ) a formidable appearance. Also (|3(]) is not too familiar when n > 2. However, 
notice that the resolvent kernel of P n on § n (M n ), with P n restricted to the orthogonal complement 
of its kernel space, is always G(x,y) = — ln\x — y\, with x,y £ S n (x,y G M. n ), and with | . | 
the chordal distance on § n (Euclidean distance in R n ), cf.|12|. More precisely, —P n ln\x — y\ = 
(l/2)(n - l)\\S n \(5 y (x) - IS™]- 1 ) on S n . Hence, all the equations @, as well as @), have a 



statistical mechanics interpretation whenever Q has a well defined sign. 

In the following we discuss the prescribed Paneitz curvature problems on S n , (^), (29), using 
the mean-field limit of the regularized microcanonical ensemble of J41[| . Nirenberg's problem on § 2 
is contained in the analysis as special case n = 2. We also remark on the equations in M n , (^0|), 
using just the canonical ensemble. 

Let us begin with the simpler (j30|). A brief moment of reflection reveals that Theorem 1, and 
its proof, have an immediate generalization to the Paneitz equation (|30|). In (|l]) we then have to 
set G(x,y) = — ln\x — y\, x,y £ M n . We also replace \K\ in by and 2h by n/i, where 

instead of a harmonic function as in Theorem 1 now h is a non-constant element of the kernel 
space of (-A) n / 2 on R n , to which we may want to refer as higher harmonic function. Moreover, 
Q shall have well defined sign and satisfy analogous integrability conditions as ( |l4| ) and ( |i~5| ) in 
Theorem 1. Finally, the critical /3 = — 4 for the canonical ensemble changes to (3 = —2n, and the 
corresponding critical integral Gauss curvature k = 4ir in Theorem 1 changes to a critical integral 
Paneitz curvature q = (n — l)!|S n |, where q = J Rn Q(x)e nu<yX ^dx. These numbers are also the sharp 
constants in the Trudinger-Moser type inequalities on A C W 1 and S n , [37, ||, 12]. 
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We now come to the problem on S n . Since S n is a compact manifold without boundary we do 
not need any delicate unbounded-domain estimates of the sort needed in [^] or [42]. The technique 
of [41] carries over to S n without major changes. In ([!]) we simply have to set G(x,y) = — ln\x — y\, 
x,y G S™. We also let F^ N \x) = Nf(x), with / some continuous function on S n satisfying 
Jsn f(x)dx = 0. Notice that not all / will correspond to a Paneitz curvature function. 



Our regularized microcanonical probability measure on (S n ) N is of the form El 



(N,e,a) 



(dxi...dx]y) = — exp 



-N— (s-— ijWV 
2a 2 V N 2 J 



dx\...dxN , 



where a > and e are fixed real numbers, and 



Z(N,s,a) 



exp 



(S n y 



1 



1 



-N—(e-—HW 
2a 2 V N 2 



dx\...dxN 



(31) 



(32) 



The Hamiltonian is given in (||), with the identifications of G and F mentioned above. The 
microcanonical ensemble at fixed N is obtained in the limit a — > in (|3l|), giving a delta measure 
concentrated on = E}, with E = N 2 e, as can be easily verified using geometric measure 

theory. 

Let P((S n ) ) denote the probability measures on (S n ) . For any G N, we define the entropy 
of qn S P((E> n ) N ) relative to the normalized uniform measure \§> n \~ N dxi...dxN by 



S(g N ) 



(S n ) N 



PArln(|S n | prA dx\...dxjsi , 



(33) 



if qn is absolutely continuous w.r.t. uniform measure on (§") , having density pn, and provided 
the integral on the r.h.s. of ( |33| ) exists; S(qn) = -co in all other cases. For Qi = g £ P(S n ), we 
define a one-particle penalized entropy functional by 

ReAe) = S(e) - ^ ( e " \J f sn y)Q(dx)g(dy) - J f(x)g(dx)j (34) 

for those g[dx) = pdx, with p a probability density, for which S(g) > — oo. We set R £jCr (^) = — oo 
in all other cases. Here, S(g) = S(qi) is the one-particle entropy as defined in (p3|). We write M £iCT 
for the set {g £ ,a} C P(§ n ) of maximizers of R £ , CT . 

By Q = (S n ) N we denote the § n - valued infinite exchangeable sequences, by P sym (p) the 
permutation-invariant probability measures on f2. According to the theorem of de Finetti - Dynkin 



1 32] , |3Q] every p G P sym (yi) is a unique convex linear superposition of product measures, see also 
[28] and [35]. By a theorem of Hewitt and Savage 1 35] , the product states £>® N are also the extreme 
points of the convex set P sym (Q). Hence, we have the extremal decomposition 



p k (dxi...dx k ) 



P(S») 



v{p\dg) g® k {dx\...dxk) 



(35) 



with pk{dx\...dxk) G -P((S n ) ) the fc-th marginal measure of p. 

As in [^TJ, we first take the limit A^ — > oo of (31) for fixed e and <r. This gives us mean-field 
theory with a two-parameter penalized entropy principle, as follows. 
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Theorem 2: For each e£l and a > fixed, one has 

lim ±-ln[\S n \- N Z(N,e,a)] = R £ , CT (f? E|(T ) (36) 

iV — >oo IV 



urai/i £ £)(T G M £i(7 . Moreover, (31) has at least one limit point in the corresponding subset of P sym (VL), 
convergence understood for all the marginals in the sense of Kolmogorov, J^, here weakly in LP , 
p < oo. The decomposition measure v(p^ £,(7 ^\dg) of any limit point p,( £ > a > is concentrated on M £)Cr . 

The subsequent limit a — > now gives the anticipated mean-field variational principle for the 
microcanonical entropy. We denote by L}' + (S n ) the subset of the positive cone of L 1 (S n ) whose 
elements p satisfy J* § „ pdx = 1, and by l}^(S") the subset of L}' + (S n ) for which 

\ I [ G(x,y)g(dx)g(dy) + [ f(x)g(dx) = e, (37) 

with g(dx) = p{x)dx. Let £o(/) denote the minimum over P(E> n ) of the functional on the left side 
of ©. 

Theorem 3: Fix eel. 

Part 1. Let e > Eq. Then the limit 

s(e) = lim R £ ,a(0e,a) (38) 

cr— >0 

exists and satisfies the variational principle 

s(e) = max js(£>)| g(dx) = pdx; p G ^i^j • (39) 
Lf e > Eq, all maximizers p £ for ij^t) satisfy the Euler- Lagrange equation 



p £ {x) = exp [j3 



Pch - / G(x, y)p e {y)dy - f(x) 



(40) 



where (3 and p c h are real Lagrange parameters for the constraints p G ^i-t ■ F° r £ = £ o> 
maximizer(s) solve the free boundary problem obtained from in the limit (3 — > +oo. 

Furthermore, let M £ denote the set of maximizers p £ for (jg^). Let p^ be a weak limit point, 
as o~ — ► ; of the measure pj- £,cr > . Then p^ E P sym ({l), and its decomposition measure v(p^\dg) 
is concentrated on M e . 

Part 2. Let e < Eq. Ln this case lim^^o ^e,a(Qe,cr) = — oo. 

The proofs of Theorems 2 and 3 are nearly verbatim copies of the proofs for two-dimensional 



domains given in 41], with a few functional analytical differences. Details will appear elsewhere. 

Our ( |4C)|) is the dual equation to (|28| ) on S n . The differential geometric problem is to find such / 
for which ( [Io|) has a solution with (3 = — 2n, in which case (n — 1)!|S™| exp (2nf(x)) can be identified 
with a prescribed Paneitz curvature functions Q(x). Clearly, only such Q with sign (Q) = +1 can 
be found this way. The parameter (3p c h, which accounts for the requirement that p is a probability 
density, is simply absorbed in u(x), but the parameter (3 in ([4C|) is implicitly determined by the 
choice of e. Hence, the problem becomes: "Find all / for which the map e i— ► [3(e) takes the value 
-2n." 
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One obvious such / is / = 0, in which case the equation ( [40]) with j3 = —In is completely 
integrable [12, Namely, (40) is then invariant under the full conformal group on S n , having a 
unique solution — up to rotations and reflections on S n — for each e > Eq. 

The problem to find admissible / ^ may not appear any simpler than the original Nirenberg 
problem and its generalization to higher n, but now thermodynamics comes to the rescue. It is 
straightforward to show that s(e) is continuous and piecewise differentiable and that Ran(<9 e s) is 
connected. As in ordinary thermodynamics, so also here we have the identification 

= d £ s(e) (41) 

wherever the derivative is defined. We have d £ s(e) — ► +oo as e — > Eq, except when / = 0. Moreover, 
for all / ^ we have s(£oo) = and <9 e s(£oo) = 0, where £«, = (1/2) L n J* §n G(x,y)\S n \~ 2 dxdy. For 
£ > £ooj we have d £ s(e) < 0. Therefore, it suffices to solve the simpler problem of finding those 
/ for which d £ s(e) < —2n for some large enough e > £oo- Using physical intuition as guidance, a 
little reflection reveals that a solution to the generalized Nirenberg problem should exist whenever 
/ allows particles to cluster in at least two spatially separated regions when (3 < 0, or, in technical 
language: 

Conjecture 4: Any f S C°(S n ) satisfying J gn f(x)dx = which has at least two isolated maxima 
can be identified with a Paneitz curvature function by Q(x) = (n — 1)!|S"| exp (2nf(x)). 

A variant of Conjecture 4 for §> 2 with somewhat stronger regularity assumptions on the Gauss 
curvature has been proven with PDE techniques in see their Theorems I and II. Our micro- 
canonical mean-field limit now offers a new perspective on proving a result like Conjecture 4 - 
simultaneously in all dimensions. I hope to report on the details of such an effort in some future 
publication. 

As a final remark, I mention that also the ground state problem e = eq, while only indirectly 
relevant to our conformal geometric problems on S n , is of quite some interest, in other contexts; 
sec [ 



6J. 
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